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Aim of this study is to show a way to unifiedly understand the origin of repulsive core of 
baryon-baryon interaction and the origin of universal repulsion of three-baryon interaction 
whose importance has recently become noticeable. "Universal" means its effect independent 
of flavor and spin of baryons. As shown in recent lattice QCD calculations, the confinement 
mechanism in the single baryon is realized in the Y-shaped structure of the strings em- 
bodying the squeezed color flux-tubes that consists of a junction at the central position and 
three oriented strings linking the quarks with the junction. Here we adopt a string-junction 
model suitable for description of such baryonic structure, in order to find out substantial 
features that the confinement mechanism persists in multi-baryon systems when baryons 
closely approach to each other. First we clarify an originating mechanism of the repulsive 
core of baryon-baryon interactin, and next by extending this approach to the three-baryon 
system, we derive a universal short-range three-body repulsion. Key concept in the study 
lies in the recognition that such short-range repulsion appears as latent effect implying the 
energy necessary for full overlapp of baryons, that is, the energy for creation of the junction 
and anti-junction pair which leads to specific exotic baryon states of gluonic excitation. The 
short-range interaction derived in such a way is universal, because the effect comes from 
the color degrees of freedom. This three-body repulsion prevents the equation of state from 
the dramatic softening due to hyperon-mixing in neutron-star matter, and thus we can dis- 
olve the difficulty that otherwise theory contradicts with the observations on the masses of 
neutron stars. 

§1. Introduction 

In recent studies of neutron star (NS) matter with a hyperon-mixed core at 
high density, Takatsuka et al.^^ have shown that universal repulsion of three-baryon 
interaction (3BI) is indispensable to avoid the dramatic softening of equation of 
state (EOS) due to the hyperon-mixing, which leads to the contradiction to the ob- 
servations on the masses of NSs.***-* Here "universal" means that the repulsion is 
independent of flavor and spin. This recognition is significant, although phenomeno- 
logical, because the statement holds almost independently of the models describing 
NS matter. Importance of the recognition is similar to the case of the repulsive core 
model for nuclear force proposed by Jastrow in 1951,^^ who succeeded in explaining 
the high energy nucleon-nucleon scattering data at that time without violation of 
the charge-independence of nuclear force. Therefore it is a very interesting problem 
to study the origin of such a universal 3BI repulsion. In three-body force, there is 



*' Preliminarily reported in Ref. f). 
**' E-mail: tama-ktn@nike.eonet.ne.jp 

*•*) ^jjg following abbreviations; NS for neutron star, 3BI for three-baryon interaction, 

EOS for equation of state, and SJM for string-junction model. 
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3BI of meson-exchange nature acting in the intermediate region as weh as 3BI acting 
at small inter-baryon distances. In this paper we concentrate on 3BI of short-range 
nature. 

Recently we have studied the problem whether or not the 3BI due to the two- 
pion exchange via /^-excitation gives rise a universally repulsive effect in high-density 
neutron star matter, and obtained the following results.^) Although this 3BI acts 
repulsively in pure neutron matter, it does not provide a repulsive effect at high 
density strong enough to prevent the EOS of NSs from the dramatinc softening due 
to the hyperon-mixing, because it does not acts on the admixed yl-particle due to its 
spin-isospin dependence. Even though there may be other possibilities of mesonic 
origin, to get universally rcplusive effects in 3BI is quite open. 

Noticing the property that the repulsive 3BI needed phenomcnologically is com- 
mon for all the baryons, we study the problem from the viewpoint of the quark 
confinement, the nonpertubative effect inherent in the color degrees of freedom in 
QCD. In treating the baryon system, we should adopt a model which well describes 
the confinement mechanism in the single baryon. In recent studies with use of the 
lattice QCD calculations, Suganuma et al.^^ have shown that the confinement mech- 
anism in the single baryon is represented as the Y-shaped structure consisting of a 
junction at the central point and three oriented strings embodying the squeezed color 
flux-tubes. The junction plays a role to neutralize the color and the strings stem 
from the junction and end at the quarks. Such a Y-shaped string-junction picture 
of the single baryon was firstly pointed out by Nambu in 1973,^^ and adopted as 
one of key structural elements in the string-juntion model developed in the lOTO's."^) 
The recent lattice QCD calculations have given the basis that this picture for the 
confinement in the single baryon is real. In the present paper we study the problem 
by adopting the string-junction model we abbreviate as SJM. 

At first, developing a previous study based on the SJM by the author in 1982,®-* 
we show a way to understand the originating mechanism of the repulsive core in 
the two-baryon system, where we assume the existence of exotic dibaryon states 
predicted by the SJM. Extending such a line of approach to three baryon system, we 
derive the universal repulsion of 3BI with short-range nature in the SJM. In the next 
step we study the effects from this part of 3BI on the EOS of neutron matter, by 
using the framework constructing an effective density-dependent two-body potential 
from 3BI. Because the resulting 3BI is universal, its effect in the NS matter does 
not change with the hyperon-mixing. Finally we show that the EOS obtained by 
including the effects from this universal repulsion of 3BI is similar to the sufficiently 
stiff EOS given by Takatsuka et al.^^ Thus we can avoid the drasmatic softening of 
EOS due to the hyperon-mixing which otherwise leads to the contradiction with the 
observations on NS masses. 

This paper is organized as follows. After a brief description on the string-junction 
model in §2, we show the way to understand the origin of the BB repulsive core in 
§3. In §4, the derivation of the universal repulsion of 3BI in the SJM is described, 
the density-dependent effective two-body potential is derived from this 3BI repulsion 
and results of numerical calculations are shown. Then in §5, the effects on the EOS 
of NS matter are numerically calculated and the results are shown. Finally several 
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related points are discussed in §6. Concluding remarks are given in the last section. 

§2. Brief description of the string-junction model (SJM) 

In this section, we briefly explain the SJM and give mass formula of hadrons 
which are used in later sections. 

2.1. Description of ordinary hadrons in the SJM 

The SJM was developed in the later half of the 1970's, following the studies of 
hadron reactions with use of quark lines and their rearrangement diagrams.^)' A 
scries of studies on the SJM by Imachi, Otsuki and Toyoda-'^^)"-'^'^) arc to be especially 
noted. Study of the present paper is performed on the basis of these works. In the 
SJM the color flux-tubes embodying the conflnement mechanism are illustrated by 
strings with orientation, as in Fig. 1. For mesons, we define the orientation of strings 
as stemming from the anti-quark (q) and ending at the quark (q). (The reverse choice 
is equally possible, as this is matter of definition.) The confinement mechanism is 
represented by formation of a closed net of oriented string(s). An ordinary meson is 
drawn by a single string as shown for M in Fig. 1. An ordinary baryon (anti-baryon) 
is drwan by the Y-shaped net with three strings stemming from a junction J (ending 
at an anti-junction J), as shown for B (B) in Fig. 1. The arrow attached to the string 
represents its orientation that means the predisposition of the order in appearance 
of q (denoted by o) and q (denoted by •) when we cut the string. This holds also 
for the case that q and/or q do not appear at the end. Following our definition, q 
appears ahead of the arrow and q appears behind the arrow. We simply write the 
ordinary meson of {qq) as M, the ordinary baryon of (qqqJ) as B and the ordinary 
anti-baryon (qqqJ) as B, where (■ ■ • ) denotes the color singlet object. 



I 1 

t fm 




Fig. 1. Pictorial view of the ordinary hadrons; M, B and B. The figures do not represent individual 
states of hadrons, but rather imply the states from which actual states are generated as Regge 
trajectories after the angular momentum projection, like the so-called intrinsic states in deformed 
nuclei. When wc superpose various rotated configulatioris in the body-fixed frame, there appear 
the states spreading smoothly. The dotted curves illustrate the region supposed in the bag 
model. 

The string structure for M has been noticed as a suitable represention for the 
color flux-tube since the early stage of QCD study. For the structure of B, Suganuma 
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and others have shown, in their recent lattice QCD calculations, that the junction 
structure of the Y-shaped strings is real.^^'^^^ This recognition has given the sound 
basis to such a view adopted in the SJM studies in the 1970's. The results of these 
QCD simulations are summarized in the energy expression of a potential form:^) 

Vgg{r) = a,gr-^ + C,g, (2-1) 

Vsg = CrsgLmin — A^g ^ -j — + C^q, (2-2) 

(73q ~ 0.89 GeV/fm ~ (Tgq, A^g ^Agg. (2-3) 

Here r for M is the distance between q and q, and Lmin for B is the minimum of the 
total lentgh of the color flux-tubes linking three quarks with J. The 1st term of V^g 
is the energy of the strings with the string tension almost equal to that of the q-q 
case; this means the universality of the string tension. The 2nd term represents the 
color Coulomb-type attraction of the onc-gluon exchange nature, as inferred from 
the relation A^g ~ ^Agg. Such presentation has been shown to be valid over the 
actual range of Lmin ~ (0.5 — 1.5) fm. The lattice QCD calculations have shown the 
profile of this Y-shaped string-junction in the action density, and the actual length 
of the strings from the junction to the quarks, Ig ~ 0.5 fm,^^) that is, -Lmin — 1-5 fm 
for the equilateral configuration of the quarks. 

2.2. Mass form,ula for hadrons in the SJM 

Usually the hadrons other than M, B and B are called exotic hadrons. This 
terminology is somewhat ambiguous, because there possibly exist many kinds of 
hadronic states, which are not simply described in the form of (qq), (qqqJ) and 
(qqqJ) but not necessarily exotic, e.g., molecular states of the ordinary hadrons. In 
the SJM, we can descriminate exotic hadrons from nonexotic hadrons; exotic hadrons 
are characterized by the appearance of the interjunction string (abbreviated as IJ) 
which connects J and J. Hadrons are constructed by three kinds of building blocks, 
that is, the quark q (anti-quark q) with one oriented string, the junction J (anti- 
junction J) with three strings and the string of the inter-junction IJ. We define 
the total number of quarks (the sum of the numbers of q and q) denoted as A^^, the 
total number of junctions (the sum of the numbers of J and J) by Nj and the total 
number of interj unctions by Njj. Then we denote the energy of each building block 
by nig, mj and m/j, respectively, where we take the same mass for q and q and the 
same energy for J and J. 

The mass formula in the SJM are given in the simplest form, as follows: 

m{Ng, Nj, Nij) = mgNg + mjNj + mijNu = uieNj - 5Nij, (2-4) 
where we use the relations. 



Nq = 3Nj-2Nij (Njj^O), mB = 3mg + mj, d = 2mq-mij. (2-5) 
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The mass niq represents the sum of the current quark mass and the energy of an 
attached string. For the hght quarks {u, d and s) the string energy is the main part 
of rUq. In this paper we treat such a case. The mass ms represents the average of 
the ground state of the single baryon. If \5\ is much smaUer than the baryon mass 
~ 1 GeV, the energies of hadrons are mainly given by msNj and spread by the 
secondary contribution —5 Njj. Estimate of magnitude of S is discussed later. For 
the time being, we proceed under the assumption of the property \d\ <C tub, which 
is proved in §6. 

The most importnat point in this mass formula lies in the property that the order 
of a hadron mass is determined by the number of the junctions as m{Nq, Nj, Nij) ~ 
Nj GeV, which we preferentially use in the following arguments. 




Fig. 2. Possible features of strings in approaching two baryons; (a) two adjacent strings change 
their connections with the nearest junctions to make the total length of strings minimum (flip- 
flop prpcess) and (b) at closer approach, two B suffer strong deformation but the strings still 
keep the Y-shaped structure.®®® 



2.3. Excitation of ordinary hadrons 

A typical order of magnitude for increase of energy of M and B is considered as 
that of the masses of vector mesons such as p and to. Meson production as fission of 
the string takes place for the elongation of a string with the order of 0.5 — 0.7 fm in 
M and B. Excitation with such order of magnitude is also possible by deformation 
of the flux-tube and by swing in a relative angle between strings, which are regarded 
as vibrational modes in a broad sense. These are the excitations keeping the basic 
structure of M and B, and the substatial features of the confinement do not change. 

At close approach of two baryons, when two Y-shaped strings begin to overlap, 
the rearrangement between two adjacnet strings takes place to make deviation of the 
total length of the strings minimum, as shown in Fig. 2(a), which is called the "flip- 
flop" process of strings.^''-* At closer approach, strong deformation of the Y-shaped 
strings necessarily occurs so as to keep the Y-structure, as illustralted in Fig. 2(b). 
Resulting excitation energy is supposed to be considerably large. 
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§3. Origin of repulsive core of baryon-baryon interaction described in 

the String-Junction model 

At present it is an open problem whether the repulsive core*) exists commonly 
in all the BB states or changes its aspects depending on flavor and spin. At least 
it certainly exists in the ^Sq and ^Si states in the nucleon-nucleon (A^A^) system. 
Although there is no sound evidence for its existence in the A^A^ P-waves and higher 
partial waves, the assumption of its existence in all the A^A^ states does not contra- 
dict with the experimental data. As for the systems involving hyperons (Y), there 
is no experimental information on specific properties of short-range NY and YY 
interactions. Most of theoretical studies are usually done following the NN case 
under the assumption of approximate falvor symmetry. At present the assumption 
that the repulsive core exists in all the states of the BB system does not contradict 
with the observations. 

From the viewpoint of the SJM, because we consider the repulsive core as an ef- 
fect originating from the confinement mechanism, the universal repulsive core (acting 
commonly for all the baryon pairs) is a natural consequence. 

3.1. Exotic dibaryon states and junction-pair states 

When two B approach very closely, the two Y-shaped strings inevitablly over- 
lap. At full overlap two B lose their identity because they cannot keep the Y-shaped 
structure and finally fuse. How does the confinement mechanism work at such situa- 
tion? To this problem, the SJM provides us with the viewpoint that the strings and 
junctions make the closed net of the color flux-tubes with six quarks at the ends by 
producing the exotic dibaryon state, which is illustrated in the (left) of Fig. 3. As it 
has A^q = 6 and Nj = 4, we denote it by L>|. ***) 

In this exotic dibaryon state Dq, there appear three J and one J, and the 
three interj unctions connecting them, which are absent in the ordinary hadrons, 
manifest the exotic character. Because the reality of J and consequently of J becomes 
firm, the appearance of the oriented interjunction connecting them (indicated by the 
arrow emerging from J and ending J) is naturally expected. Thus we consider the 
existence of quite reasonable. Although the auguments in the following done on 
the assumption of the existence of L>| predicted by the SJM are speculative, we can 
obtain meaningful consequence. In the last two sections, after presentation of the 
results, we make comments on such problematic points that positive observational 
information on the existence of Dq has not been reported yet and no lattice QCD 
simulation about it has been performed yet. 



*' There are several ways to describe the short-range repulsive effect in the BB interaction. We 
simply call this effect by "repusive core" in this paper. 

If the so-called H-particle, a dibaryon state of the flavor singlet, would be observed as a 
spatially compact bound state, this assumption does not hold. But at present the H-particle has 
not been observed. 
***'> We denote a hadron H with Ng and Nj by H^'' . 
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Fig. 3. (LeftjString-junction structure of exotic dibaryon Dq. (Rightj Formation of BBBB through 
fission of three interjunctions. 

The mass of Dq is expressed with use of Eqs. (2-4) and (2-5): 



The energy of Dq is higher than that of the ordinary state of 2B {2mB) by £. 
Dq can not turn back directly to BB. In order to turn into a system of the ordinary 
hadrons, the cut of three interjunctions is required. In doing so, as shown in Fig. 3 
(right), at the cut places the quark pairs qq are created and then the baryon-pair BB 
is created. As a result the system turns to BBBB. Its energy is higher than 2mB 
by about 2ms, and therefore we can take £^ ~ 2 GeV. This is simply understood 
from the fact that Nj increases by 2, compared with the ordinary 2B, according to 
the mass formula. 

In the description mentioned above, the exotic dibaryon Dq is considered as to 
be formed through the virtual creation of the BB pair. The exotic dibaryon Dq is 
also formed through another process. Such an example is shown in Fig. 4. In the 
(left) of this figure, situation proir formation is illustrated: In order for two parallel 
strings with the same orientation belonging to different baryons to overlap, one of 
them is excited by formation of the string loop accompanying the junction pair. In 
the (right) of the figure, situation just after formation of Dq is illustrated: Two 
strings with opposite orientations transform into two hair-pin like strings. The 
Dq is surely formed, although some of the strings are strongly deformed. In this case 
too, since the junction number Nj increases by 2 for the formation of one junction 
pair, the excitaion energy is about 2 GeV. Such a loop-type excitation of one string 
has been studied in a recent work on the pentaquark problem by Suganuma and 
others. ^^^'^'^^ In lattice QCD calculations in a study of the single baryon, gluonic 
excitation has been shown to appear at about 1 GeV and to be regarded as a global 

*^ In the isolated BB, mcson-cxchange effect works attractively mainly duo to the ui and a, 
but here such effect is not considered because it is uncertain in many-baryon configuration under 
consideration. 



m{Dl) 



Grug + 3m J + mj + 3mij = 2(3mq -|- mj) + mj + mj + Smij 
2mB + £, 
2m J + Smij. 



£ 



(3-1) 
(3-2) 



8 



R. Tamagaki 




Fig. 4. Illustration of an example for another process to create _Dg ; (left) an excited string appears 
to creat a string loop and a J J pair and another adjacent string is in parallel with the same ori- 
entation which is opposite to the orientation of the strings in the loop and (right) rearrangement 
takes palce between two strings with opposite directions, leading to the formation of Dg. 

excitation of the whole Y-type flux-tube system. On the contrary, the excitation 
shown in Fig. 4 is that of a single flux-tube and appears in the energy near £ ^ 2 GeV 
for the BB-iype excitation, apart from the deformation energy of the string. 

The BB pair turns into a series of states commonly involving the J J pair, as 
shown in Fig. 5; M4 and Mfiexotic mesons), and states involving no quark S^ia kind 
of glue balls). We note that, in later descriptions, £ ~ 2GeV plays a role of the key 
number. The loop in Ml created by fusion of the q and q in M4 is the same with 
the one appearing in Fig. 4. So we may call the excitations mentioned above simply 
the BB excitation or the J J excitation. 

This series of states with the baryon number=0 and Nj = 2 contain commonly 
the J and J which we call the junction-pair. Here they are denoted by the symbol Ai. 
Because of irrelevance to the number of quarks, they are the states inherent in the 
gluon degrees of freedom. The energies of M. are distributed about the central value 
of about 2 GeV with a spread of several hundreds MeV, as far as we are concerned 
with u, d and s. (For details, §6.1 should be referred to.) This means that there 
exist a group of states with the baryon number=2, consisting of and [BB+ M.\, 
in the energy region higher by about 2 GeV above the ordinary BB. 



BB ^° 

Fig. 5. A series of states with the baryon number =0 formed by fusion of q and q from SB, which 
we call junction-pair states, denoted as M. 
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3.2. Repulsive core of the BB interaction 

We can paraphrase the statement mentioned above as follows. In order for the 
confinement mechanism to persist when two baryons overlap fully, it is necessary to 
form the string-junction net where the arrows attached to the oriented strings are 
closed. It requests the excitation energy of <S ~ 2 GeV, which makes the creation of 
BB pair (more generally the junction pair) possible. If the kinetic energy of BB in 
the center of mass system is much smaller than ms (e.g., in the elastic scattering 
region), the system is apt to avoid such a large excitation. The potentiality that two 
baryons reside apart to get rid of such expense is equivalent to the existense of a 
high potential barrier between two baryons, namely, the repulsive core. Its height is 
of the order of £^ ~ 2GeV, and its range is the relative distance where two baryons 
begin to overlap largely. 

The potential mentioned above may be defined in a sense of the adiabatic ap- 
proximation usually adopted to get a static potential, where the relative kinetic 
energy of two baryons is taken into account in solving the Shrodinger equation. In 
principle, we should treat the energy described in the quark-gluon level of QCD, but 
here we use the approximate value obtained in the SJM, abbreviated as E^^'^^\ We 
define an interaction potential dominating in the core region to be used in the baryon 
level, by the difference between the energies E^^"^^^ at two relative distances of BB, 
^12 = 1^*1 — ^"21 and ri2 — ^ oo, where r\ and r2 are the positions of two baryons: 

y(core)^^^^) = £;(^-^^)(ri2) - E^^-^^\n2 ^ oo). (3-3a) 

For the SJM, since we can give the values of E^^"^^^ only at the two limiting cases, 
at ri2 <C is aud ri2 ^ is-, compared with the string length of the ordinary baryon 
is ^ 0.5 fm, the magnitude of V^^^'^^'> at closest approach is obtained as 

y(core)(^^2 < is) ~ E^^J^\r^2 < 4) - £^^^-^^^(^2 » 4) (3-3b) 
~ E{dI) - 2mB ^ 2mB + £- 2mB ^ £■ (3-3c) 

This is regarded as the interaction of short range in the hadron level, and 

is incorpolated into the two-baryon potential, together with the meson-exchange 
potentials of the long and intermediate range. 

The view that the core-like short-range repulsion is a latent effect reflecting the 
energy needed for the overlapp is quite natural because it is familar in the atomic 
and nuclear systems. The repulsive core between two atoms (typically between 
two closed-shell atoms), which is often approximated by a hard-sphere potential, is 
understood as originating from the large excitation energy of electron states which 
is needed for the overlap of the electron clouds. This repulsive core is "strong"; it 
is so steep that its core height becomes by several order of magnitude larger than 
the strength of attractive pocket just ouside the core, if a soft core potential is 
used. The repulsive core used in the effective potential between two nuclei (typically 
between two light closed-shell nuclei, e.g., between two a-particles) is understood as 
originating from the excitation energy of nucleon states required for overlap of two 
nuclei. This repulsive core is "moderately strong", as its height of ~ 150 MeV is 
larger by one order of magnitude than the depth of attractive pocket (~ 10 MeV) 
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for the a — a case, although a hard core potential may be used phenomenologically 
for low energy scattering. *^ The BB repulsive core of the ~2 GeV height with the 
~ 0.1 GeV attractive pocket (in the iS-wave) belongs to the "moderately strong" 
class. 

These features of the repusive core in the atomic and nuclear systems are at- 
tibuted to the effects of the Pauli-exclusion principle. The related internal degrees 
of freedom is the spin of the electron in the atomic system and the spin-isospin of 
the nucleon in the nuclear system. As for the repulsive core of the BB system, the 
related internal degrees of freedom is the color. Although we see similarity among 
three cases, it is to be noted that the confinement mechanism is essential in the BB 
case, and the exclusion principle relevant to the color degrees of freedom alone is not 
essential because the quarks have many internal degrees of freedom. 

The range of the repusive core originating from the confinement mentioned above 
is smaller than the relative distance where two Y-shaped structures begin to touch. 
If the string length of the Y-shape in the ground state of B is taken as ig ~ 0.5 fm 
correspoding to -Lmin — 1-5 fm in Eq. (2-2), the half overlap of two B corresponds to 
the relative distance ~0.5 fm. Thus the range of about 0.5 fm is a reasonable choice 
for the repulsive core derived in the SJM. In a previous work in 1964,^^^ where 
we considered "structural core" of the NN interaction from the viewpoint of the 
composite model of hadrons, we have chosen the shape of the A'^A'' repulsive core as 
the gaussian form, implying that the repulsion grows as the overlap of two nucleons 
becomes large. On the same standpoint, the present author constructed a realistic 
A^A^ potential with the gaussian repulsive core, called OPEG.^^^ In a representative 
version of OPEG potential, OPEG-A, the strength and the range of the repulsive 
core are taken as 

Vc{r) = V^ex.-p[-{r/r]c)^], with = 2 GeV, rye = (0.45 - 0.5) fm. (3-4) 

The repulsive core V^^°^^^ obtained in the SJM has the feature well corresponding 
to Vc of OPEG-A. Thus we can understand the origin of the BB repulsive core in 
the SJM, as the manifestation of the confinement mechanism in the BB system. 

In the notion of the state vector, the present context is interpretated as follows. 
The state vector of the system l^') is a superposition of the part spanned by the 
states of two ordinary baryons, denoted as Pl'l^), and the part spanned by the Dq 
plus [BB + Ai] states, denoted as Q\^), where P and Q are the projection operators 
{P+Q = 1). The latter states Q\^), whose existence is limited to the core region and 
at high excitation, has the latent influence on the former state P\^)- Such influence 
appears as suppression of the BB relative wavefunction at very small distances, 
corresponding to the effects by the repulsive core. 

It is not simple to show logical connection between the core-like repulsion shown 
in the adabatic approximation and the effective BB interaction described in the P\^) 

This repulsive soft core in the potential form is understood as a reflection of the forbidden 

relative states which should vanish under the Pauh-exclusion principle. Such a feature is described 
by using the orthogonality condition to the forbidden relative states or by using a strong nonlocal 
potential. 
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influenced by the existence of Q\^). First of all we notice that substantial effects 
appear in two ways. One is the effect of the orthogonality between P\^) being in 
the scattering state and Q\'l^) with spatially compact sructurc and high excitation 
energy. The other is the channel coupling effect. The resulting core-like repulsion 
in the adiabatic approximation means that the former overwhelms the latter. The 
reasoning is as follows. 

1. Effect of the orthogonality 

In the scattering problem, as the energy is fixed to the incident energy, the 
observable effect appears in scattering phase shifts in the asymptotic region. 
In the problem under consideration, the influence from Q\^) is rather special, 
because its existence is limited in the inner-most region of the BB relative 
distance (core region), r <h 0.5 fm. We notice that the primary inflTicnce 
from the orthogonality to Q\^) brings about the nodal behavior of the BB 
relative radial wavefunction of P\^) at the core region, which is to be described 
by a strong nonlocal interaction. The reason is that such radial nodes do not 
mean attraction to produce bound states but the almost energy-independent 
nodal behavior gives repulsive effects, because the amplitude of the oscillatory 
wavefunction is suppressed and the outermost node plays a role equivalent to the 
boundary condition given by a core- like repulsion with the radius at r ~ 6, in off- 
resonance situation.*) If there exist many such states of Q\^), the wavefunctions 
are strongly sTipprcssed in the whole core region and the situation becomes as 
if the hard core exists. If some attractive potential works outside the core 
region, even though the scattering phase shifts show attractive feature at low 
energies, repulsive feature appears at high energies because repulsive-core effects 
are strong enough to overwhelm the attractive effects. The repulsive core is to 
be considered as a simple local potential form in such a limited sense that it 
represents only the suppressed feature of the wavefunction given by the nonlocal 
interaction in the core region mentioned above. **^ 

2. Channel coupling effect 

In addition there appears a dynamical effect due to the channel coupling through 
the interaction matrix elements (i.m.e.) between P\^) and Q\^). Generally the 
channel coupling effect is considered to give attractive effect on the channel with 
a lower energy, in ordinary situation of bound states. Since the states oi P\W) 
and Q\^) under consideration have substantially different structures with the 
large energy difference, the i.m.e. are not large, and the channel coupling effect 
seems moderate. We cannot enter further into this problem, because of lack of 
concrete knowledge about the wavefunctions of Q\^) in the SJM. 
We can summarize the results in the SJM approach to the BB short-range inter- 
action, by saying that the repusive core preventing two baryons from close approach 

*•* Previously wc showed such an example in the NN ^So scattering. ^^-^ 

**' We have encountered such problem concerning interaction between nuclei: The strong non- 
local interaction giving rise to the inside oscillatory wavefunction appears when we treat the inter- 
nucleus coordinate as a dynamical variable with use of the resonating group method, while the 
repulsive core with literal meaning appears when we treat the inter-nucleus distance as a variatianal 
parameter with use of the generator coordinate method. ^^^'^'^^ 
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appears as the latent effect on the two-baryon states, corresponding to the high 
excitation energy needed for the system to fully overlap. 

§4. Short-range three-body repulsion in the BBB system in the 

string-junction model 

In this section, with use of the S JM, at first we seek main aspects of three-baryon 
interaction (3BI) originating from the confinement, next give its potential form and 
strength and then derive a density-dependent effective two-body potential which is 
used to obtain its contributions to the energy of baryonic matter. Before entering 
individual steps, we briefly describe the treatment of 3BI adopted here. 

4.1. Outline of treatment of three-baryon interaction in baryonic matter 

The Hamiltonian, which enables us to understand properties of nuclei and bary- 
onic matter, is expressed in the nonrelativistic treatment as follows: 

H = J2Ki + J2v{i,j)+ J2 V^'^\i,j,k) + --- , (4-1) 

i i<j i<j<k 

where Ki is the kinetic energy of the i-th baryon Bi (including the rest mass dif- 
ference from the nucleon mass), V two-baryon interaction and thrcc-baryon 
interaction. The ellipsis means four-body (and higher) interactions, which are ne- 
glected. We take the following form for V^^^\i,j, k) consisting of the cyclic sum: 

F(3S)(^,i,fc) = ^ W{i,j;k) = W{i,j;k) + Wij,k;i) + W{k,i;j), (4-2) 

cyclic 

where W{i,j;k) represents the 3BI which arises additionally from the approach of 
the k-th. baryon to the pair of and Bj. Such a form is familiar, for example, 
in the study of the three-nucleon interaction from the two-pion exchange process 
via Zl-excitaion of the 3rd nucleon N3 intervening between two nucleons (A^i and 
N2), initiated by Fujita and Miyazawa.^^^ (We take i = 1, j = 2, A; = 3 as 
the representative one without loss of generality.) Also the repulsive three-nucleon 
interaction introduced phenomenologically by Phandharipandc and others has such a 
form.^^^'^'^) As shown below, this form is suitable for the 3BI derived in the SJM.*^ A 
primary problem is how to describe a representative one of the SJM type, W{1, 2; 3). 

Once 14^(1, 2; 3) is obtained in the SJM, next we derive the effective two-body 
potential by taking the sum of the diagonal matrix elements over the the occupied 
states of the 3rd baryon B3 in the baryonic matter. This effective two-body potential 
importantly depends on the baryon density p, and we denote it by U{1, 2; p). At this 

*^ The cyclic sum is naturally understood by viewing the situation shown in Fig. 6 of §4.2. 
The 3BI under consideration in the SJM is regarded as arising through the virtual creation of one 
BB pair in B3 approaching to the two baryons B1B2 near by. The formal structure in this case 
is analogous to the case of Fujita-Miyazawa type three-nucleon interaction, if the Zi-excitaion in 
the latter is supposed to be replaced by the one BB pair excitaion in the former. The expression 
of Eq. (4-2) implies a wider situation and is applicable to other cases not restricted to the cases 
mentioned above. 
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step we take into account the effects of the short-range correlation (abbreviated as 
SRC) coming from the interactions between Bi — B^ and between B2 — B^, while 
the SRC between Bi — B2 is included in the final step of the energy calculation. 
Then we can calculate the energy of the baryonic matter E in the framework of the 
reaction matrix theory, by adopting V{i, j) + U {i, j; p) as the two-body potential. In 
this way the symmetric sum of V^^^\i,j, k) in Eq. (4-2) is automaticUy included in 
the sum over all the two-baryon pairs. The interaction energy E'int is given by 

occupied 

£^int = 2 E (0a(l)</)/3(2)mi,2) + C/(l,2;p)|V',,^(l,2)-V'a,/3(2,l)), (4-3) 

where ^q:(1)^/3(2) is the unperturbed (plain-wave) function of the baryon pair and 
'i/'a^^(l,2) is the pair wave function in the matter which includes two-body corre- 
lations. In addition to the contribution from the usual two-body potential V , the 
effects of the 3BI in the SJM is given as the contribution from the the effective 
two-body potential U . 

4.2. Exotic tribaryon states in the BBB system and strength of three-baryon force 

Here in the SJM, we study the representative part of the 3BI, W{1, 2; 3), which 
arises additionally when the 3rd baryon B3 approaches to the adjacent two baryons 
{Bi and B2). We show that a short-range three-body repulsion originates from 
the confinement mechanism, by extending the reasoning in deriving the universal 
repulsive core of the BB interaction in the previous section to the three-baryon 
system. 

In order for baryons to overlap fully, from the viewpoint of the SJM, the strings 
and junctions representing color flux-tubes should be rearranged to hold the con- 
finement in such situation. In the BB system, it becomes possible by creating one 
BB pair. Here we consider the extension of this idea to the BBB system by taking 
two steps. (1) Wc imagine the situation that Bi and B2 are close enough to form 
the exotic dibaryon state Dq{1,2) and the third baryon B^ approaches to Dg(l,2), 
as illustrated in Fig. 6(a). (2) One BB pair is virtually created between them as 




D|(l,2) 



D^(l,2) B B 



Fig. 6. Illustration of formation of the exotic tribaryon in two steps; (ajapproach of B3 to the exotic 
dibaryon 2) formed by adjacent Bi and B2, (b) virtual creation of one BB pair between 

them, and (c) formation of the exotic tribaryon Tg after the fusion of the strings of q and q. 
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in Fig. 6(b), and after the fusion of the strings of q and q the system turns to the 
connected string-junction net with Nq = 9 and Nj = 7 with increase of three inter- 
junctions instead of three qq-pair strings, as illustrated in Fig. 6(c). This state is the 
exotic tribaryon state denoted by Tg . 

The energy of Tg is given using the mass formula as follows: 

m{Tg) = 9mq + bmj + 2m j + Qmij = 3(3mq + mj) + 2mj + 2mj + 6m/j 

= SmB + 4mj + 6mij = Sms + 2S. (4-4) 

The mass m (Tg^) is higher by 2f ~ 4 GeV than the energy of the ordinary BBB 
system, where £ is defined in Eq. (3-2). In this energy, its half {£ ^ 2 GeV) is 
attributed to the formation of Dq{1,2). Therefore the energy arising additionally 
from the approach of S3 is £^ ~ 2 GeV. This is clear, because the BB is created 
and the junction number increases by 2. Next we consider the 3BI based on the 
SJM, whose strength is this additional increase of energy (f ~ 2 GeV) at the closest 
approach of B^ to Dg(l,2). 

We write the 3BI in the SJM (abbreviated as SJM-3BI), following the reasoning 
adopted in the two-baryon case. A representative part of the SJM-3BI denoted by 
W{1, 2; 3) is given as the difference between the energies at the two situations taken 
along the "path" which leads to the formation of Tg at the close approach of B^ to 
the exotic dibaryon Dq{1,2). These energies are to be described in the quark-gluon 
level of QCD, in principle, but here we use the approximate values obtained in the 
SJM, denoted by £;('^-'^*^)(^3; D^(l, 2)), where ^3 is the relative distance between the 
position of B3 and that of the center of mass of D^. We define the magnitude of 




Fig. 7. Pictorial view of the exotic tribaryon and its preformation stage; (leftjstring-junction net 
of the tribaryon TJ, {rightjE BB{BB){BB) states arising from the fission of the interjunction 
strings, and (bottom) an example illustrating one of many possible configurations for full overlap 
of BBB, where the dotted area indicates the formation region of a string-junction net. Such 
view is the unfolded-sheet drawing of the tribaryon having three-dimensional spread. 

@@@ 
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W{1,2;3) as 

T^(l,2;3) = E^^'^\^s;Dl{l,2)) - e(S'^\^s ^ oo;i?|(l,2)), (4-5a) 

For the SJM, since we can give the values of E^^"^^^ only at the two limiting situa^ 

tions, namely, those at .^3 <C d-s and ^3 ^ Us (equivalent to .^3 ^ 00), compared with 
the string length of the ordinary baryon — 0.5 fm, the magnitude of W{1, 2; 3) at 
the closest approach of B3 is obtained as 

T^(l,2;3):^ E^^-'^\i,, « ls:Dt{l,2)) - E^^-^^^Xi,, » Dt{l,2)) (4-5b) 
~ E{Tl) - {ruB + E{Dl)) ~ 3ms + 2E - {tub + 2ms + £)^£, (4-5c) 

which is considered as the strength of the SJM-3BI. 

W(l,2;3) regarded as the interaction in the hadron level is a genuine 3BI in 
the following sense. In L'g(l,2), the original baryons {Bi and B2) have lost their 
identity of the baryon as shown in Fig. 6(a), although Dq{1,2) is formed from B\ and 
B2 through the virtual creation of the BB pair. Therefore, the interaction between 
-D|( 1,2) and -B3 in the configuration leading to the formation of Tg is an intrinsic 
three-body effect which is not reduced to any effect due to the two-body interaction. 

The state vector of the system \^) can be written, with use of two projection 
operators P and Q satisfying P + Q = 1, as the sum of the state vector P\^) 
which describes the states of three ordinary baryons and the state vector Q\^), 
which describes the states arising from DQ{i,j) plus Bfc and leading to Tg (denoted 
by Qt\^)) and the higher states with larger junction numbers. Restricting Q\'l^) 
to Qt\'1^), we can write Qt as Qt = Q3 + Qi + Q2, where Q3 projects out the 
states brought about along the "path" leading to Tg formed from L>g(l,2) and the 
intervener B3. Qi and Q2 are given by its cyclic permutations. This means that 
three Qk\^) spans different Fock spaces. The reason is that, for k = 3, once Dq{1, 2) 
is formed from Bi and B2, the simultaneous formation of Dq{2,3) and/or -D|(3, 1) 
does not take place, because such a double (also triple) Dq configuration cannot 
form a closed net of the oriented strings and junctions. The latent effect from the 
existence of Qsl'l^) on the three-baryon state P\^) appears as the repulsive 3BI, 
W{1, 2; 3). In the same way, combination of D^(2, 3) and Bi {Dq{3, 1) and B2) leads 
to Qi\^) iQ2\'^)), which brings about W(2,3;l) (1^(3,1; 2)) on P\>1^). Thus the 
sum, Qt\^) = iZiQkl^)^ provides the SJM-3BI as the latent effect on P\W), which 
is given by the cyclic sum of W{i,j; k) as in Eq. (4-2). 

The exotic tribaryon states can be attained also in one step from the ordinary 
BBB system by creating two BB pairs as illustrated in the (right) of Fig. 7. Fusion 
of the strings of q and q leads to the interj unctions, and the BBB{BB){BB) system 
turns into the exotic tribaryon Tg , shown in the (left) of Fig. 7, This state is the 
same with the one denoted by Ag in ref. lib). 

4.3. Effective two-body potential from three-baryon interaction in the baryonic mat- 
ter 

In order to evaluate the effect of the SJM-3BI in the baryonic matter, we employ 
a practical way, called the effective two-body potential (2BP) method, which enables 
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us to calculate the effects from 3BI in a form of the additional contribution to the 
two-baryon potential. This method was developed by Hokkaido group in the study 
on the role of the 3BI due to the two-pion exchange via Z\-excitation in nuclear 
matter, and used in recent works on 

3BI.32),33) This is a way to obtain the effects 
of 3BI approximately, in place of full scale Bethe-Faddeev calculation with 3BI. The 
procedure in calculating the effective 2BP from the SJM-3BI are as follows. 

1. First of all, effects of the 3BI appear under the influence of the primary inter- 
action, the two-baryon interaction V, that is, its main effects are obtained with 
use of the wavefunction distorted by V. In baryonic matter, the distortion from 
the plane wave is remarkable at small inter-baryon distances, and is called the 
short-range correlation (abbreviated as SRC). 

2. Applying the method, we derive an effective 2BP, [7(1,2;/)). from W'^(l,2;3) 
acting in the baryonic matter, by taking the sum of the diagonal matrix elements 
of (i.e., by taking the average) over the occupied states of the 3rd baryon B^. 

3. By the procedure of taking this average, because the states of S3 do not change, 
the effect from W{1,2;3) results in the additional effect on the interaction 
between Bi and B2. 

4. A remaining problem is to choose a functional form of W{1, 2; 3) with use of 

three spatial variables (ri, r2 and r^), which are available here because of the 
universality (independence on flavor and spin) of the SJM-3BI. 



„ (occ) 

C/(l,2;p)= / dr3^(</.^(3)|T^(ri, r2 ■,rs)\cf>^i3))firi - rs)firs - r2), 
•' 7 

= P j drzW{ri, r2;rz)p(ri-r2,)f{r2,-r2), (4-6) 

where /(r, — Vj) is a SRC function which includes the suppression of the relative 
wavefunction by the repulsive core for B^ — Bj pair. Here we assume W being local. 
4>'y{S) is the plane wave state, whose spatial part is (/>'y(3) = exp(jq^ • r^)/^/l2, Q 
being the normalization volume. The sum over the states (j)^ simply gives the baryon 
density p, because of the universality of the SJM-3BI. 

The same procedure is apphcd to W{2, 3; 1) and W{3, 1; 2) leading to the effective 
2BPs in the baryonic matter, U{2^ 3; p) and ?7(3, 1; />), respectively. 

For choice of the functional form of W{1, 2; 3), at best, we can take into account 
its property that W^(l,2;3) is the interaction mediated by B^ between Bi and B2 
and is of short-range nature. In this context we take the simplest form, the product 
of a function of Bi — B^ distance and a function of B2 — -B3 distance, as follows. 

W{ri,r2, ; ra) = Wq g{ri - r^) g{r2 - r^), (4-7) 

where g' is a function such that g{rij ^ 0) ^ 1 and vanishes at rij ^ ig and Wq ~ £. 
Here, for convenience of calculations, we choose a gaussian form for g{rij) with the 
range of 77c = 0.45 - 0.50 fm in Eq. (3.4). 

g{ri - Tj) = exp(-A4), A ~ 4" - (^-^ " ^■^) (^'S) 
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Denoting the Fourier transform of g{r) as G{q), 

g{r) = J^fdl (4-9) 

G{q) =4Tr [ drr'^e-^'' jo{qr) = (7r/A)^/2g-gV4A_ ^^.-^q^^ 



After the transform of momentum variables, [/(l,2;p) = U{ri2',p) is expressed, as 
fohows: 

Uir^2■,p) = ^Jdq e-'^-("^) J dq, h{q - q^)G{q^) J dq^ h{q 

Here h{q) is the Fourier transform of the SRC function squared /^: 

fir) = j dq hiq)e-"^-^ 

It is to be noted that U{ri2; p) is strongly density-dependent. 

Without the SRC, we have h{q) = S{q), and then the integral becomes 

Uo{ri2;p) = -^Jdq e-«--(G(Q))2 (4-13) 

= pW^o(vr/2A)3/2exp(-ArV2). (4-14) 

Hereafter we write ri2 simply as r.*^ 

Now we calculate the effective 2BP ([/). To take into account the SRC for -B1--B3 
and for -B2--B3, we calculate the two integrals with the same form in Eq. (4-11). In 
doing so we need the SRC function / extracted from the solution of the reaction 
matrix equation in nuclear matter. 

The SRC functions change state by state. At low density of neutron matter, as 
the wave contributes mainly, we can use the SRC function of this wave denoted as 
fsir). As the density goes higher, this approximation is not allowed, and we should 
add both the contributions from ^Sq wave and ^-P7=o,i,2 wave. The latter depends 
on J due to the effects from the noncentral potentials. For energy quantities, the 
contributions appear in the weighted average over ^Pj=o,i,2, which is proportional 
to (2 J + 1). We define the ^P-wave SRC function: 

/p(r) = {fir; ^Pj=o) + 3/(r; ^Pj=i) + 5/(r; 'Pj=2)}/9. (4-15) 

As the ^Sq -wave and the average P wave contribute with the statistical weight 
relevannt to the spin, we define the weigted averge of the SRC function as 

/wa(r) = {/5(r) + 3/p(r)}/4. (4-16) 

*'> Making estimate of Uo, we have Uo{r; p) = 84(p/po)e"^(''/*™^' MeV for A = 4fm-2 and Wo = 2 
GeV. This value is (5-3) times larger than that with SRC at r ~ (0 — 1) fm. So the SRC is important 
for calculating effects of three-body force. 



-92)^(92)- 
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Fig. 8. Effective two-body potential U{r;p) 
from universal repulsion of SJM-3BI. Num- 
bers attached to the curves show p/po- Pa- 
rameters of the force range A and the SRC 
function / are indicated in the figure. 
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Fig. 9. Density dependence of effects from 
SJM-3BI on the energy/baryon in neutron 
matter. The lower three curves show con- 
tributions from SJM-3BI. The upper two 
curves show the total energy/baryon in 
neutron matter including SJM-3BI, where 
two-body (V^) contribution and EOS of 
TNI6u (filled circles) are taken from ref. 
9). 



These SRC functions are constructed using the wave functions obtained for OPEG-A 
potential^^^ in the reaction matrix equation in the neutron matter. 

In numerical calculations, we obtain the effective two-body potential from the 
SJM-3BI, using /s'(r) or fwa.{r)- The damping of fs{r) at small distance is strong, 
but that of fp{r) is moderate. Also the damping of fwa.{f) is not so strong. 

We denote the Fourier transform G modified by the SRC effect by GsRc(<z): 

GsKc{q) ^ j dq, h{q - q^)G{q^) (4-17) 

poo 

= 47r/ drr^flir)e-^'''joiqr), (4-18) 
Jo 

where h is substituted by its Fourier transform h{q) = f drf1{r)e'x.p{iq ■ r)/(27r)'^. 
G{q)sRC differs from G{q) by including the SRC factor in the integral. Then 

the effective two-body potential is given as 

Uir;p) = dqq^j,{qr){G^^c{q)f- (4-19) 
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In this paper the strength of the SJM-3BI is taken as Wq = 2 GeV. Results of 
numerical calculations for U{r;p) are shown in Fig. 8 at the density p/po = 1 — 5, 
where po = 0.17 fm""' is the nuclear density. The effective two-body potential 
obtained from the SJM-3BI is a universal repulsion with the range of about 1 fm. 
Although it is of the order of 10 MeV around the nuclear density po, it grows almost 
linearly as p goes high. It is rather sensitive to the spread of the function g{r); U 
becomes stronger for a wider spread. The U{r; p) obtained from the repulsive term 
of phenomenological 3BI used by Illinois group^^^ has a longer range and a weaker 
strength than those here shown. 

Adequacy of the 3BI obtained here should be examined in the energy properties 
of nuclei. We have studied the saturation curve of the symmetric nuclear matter, 
including this SJM-3BI in addition to a realistic two-nucleon potential with a repul- 
sive core and the pionic 3BI of the Fujita-Miyazawa type, and obtained a reasonable 
saturation curve.^^^ 



§5. Effects on the equation of state of neutron matter and neuron star 

matter 



We calculate the effects from the SJM-3BI on the EOS (equation of state), 
energy/baryon, in neutron matter and then discuss the effects on the EOS in neutron 
star matter. As the SJM-3BI is universal, its first order effect in the perturbation 
(without the effects of the two-body correlation in Eq. (4-3)) is simply given as 



k r p r r 

= j dkP{k){4: J drr^U{r;p)-2 J drr'^U{r; p)jo{2kr)}, (5-1) 



where is the baryon number (here the neutron number), kp the Fermi momentum 
and the 1st (2nd) term in the curly bracket is the direct iexchange) contribution. The 
function P{k) is the probability function of the relative momentum k, normalized as 
J,'- dkP{k) = IF 



2Ae 



m = - m){k/kF) + {l/2){k/kFf}. (5-2) 



The two-body correlation between Bi and B2 including the SRC effects are taken 
into accoTint by multiplying fs{r) on the ^So wave and fp{r) on the ^P wave. Then 
the contribution to the energy/baryon from the SJM-3BI (the part of the ?7-term in 
Eq. (4-3)) is given by 

+ J drr^U{r; p){fs{r) - 3fp{r))jo{2kr)}. (5-3) 



Since the effective two-body potential U has the density dependence being almost 
linear, the growth of Eu/N is stronger than the linear density-dependence. The total 
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energy /baryon is given by the sum of the contribution from the two-body potentilal 
V and the contribution from the repulsive SJM-3BI (Eu/N), as shown in Fig. 9. 
For the former we use the result by Takatsuka ct al.'^^ The upper two curves show 
the total energy/baryon (EOS at zero temperature). The bold solid curve is for 
the case of the following parameters; the range of the function g, X = 4.5 fm~^ 
(corresponding to the spread 2± 0.47 fm), Wq = 2 GeV and the /wa for the SRC 
function. In this case we obtain almost the same EOS with TNI6u (shown by • in 
the figure) without adjustment of parameters. The predicted maximum mass of 
neutron stars is M^ax — I-SMq, which is consistent with the observations. If we use 
a little wider range parameter, A = 4.0 fm~^ (corresponding to the spread=0.5 fm), 
we can provide the stiffer EOS shown by the upper bold-dashed curve, which leads 
to Mjiiax — 2.OM0. A full accout concerning the EOS including the meson exchange 
3BI will be reported elsewhere. 

Because this SJM-3BI is independent of flavor and spin, its repulsive effect is not 
influenced by the hyperon-mixing in the neutron star matter. Thus we can get rid 
of the dramatic softening of EOS, which otherwise leads to the contradiction with 
the observations. 

§6. Discussion 

Here we discuss on several related points. 
6.1. On the parameters in the SJM 

In the arguments in the previous sections, we preferentially use the key number 
6^2 GeV, without direct use of the parameters in the SJM. Here we try to estimate 
the values of the parameters (m^, mj and mij) in the mass formula of the SJM, 
utilizing this number. 

Although niB (the masss of the single baryon B) is in considrable variation, its 
source is considered as mainly due to the difference in current masses of involved 
quarks. From the viewpoint of the SJM, since the binding mechanism is mainly 
determined by the strings and junctions, we have a ground for the statement that 
the baryon mass difference is mainly attributed to the current quark masses, if 
we disregard the effects of spin-dependent interactions such as the color magnetic 
interaction leading to the color hyperfine splitting. In the following we take the 
current quark mass~ 0, suitable for the u and d quarks. Thus we take tub — ^ 
GeV, which means Sruq + mj 2± 1 GeV in the SJM. First we prove the property, 
assumed in the previous sections, that 5 = 2mq — mij used in Eqs. (2-4) and (2-5) 
is minor compared with niB ^ 1 GeV, which holds in a way irrelevant to the total 
length of the strings Lmin. Next we discuss its implications and make the estimate 
of respective parameters by taking two typical values of L^i^. 

1. Substituting the relation nij = niB — 3mq into the expression of £ = 2m j + 
Sm/j ~ 2 GeV ~ 2mB, we have 2mq ~ mij. Then we obtain the relation 
S = 2mq — mij ~ 0, namely, 6 = O(O.lGeV). This shows that \S\ <C ttib holds, 
irrespective of L^i^. 

2. The relation 2mq ~ m/j means that the energy for the fusion of q and q 
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accompanying strings is almost equal to the energy of the inter-junction string 
mij. Owing to this property, we can show that all the states in the M series 
shown in Fig. 4 have the energy of about 2 GeV. Accordingly they can transform 
to each other without large energy transfer. 
3. By taking two typical values of Lmin ~1.0 fm and 1.5 fm, we give some estimate 
for the SJM parameters. Lmin —1.0 fm (1.5 fm) corresponds to the string length 
from the junction to the quark, ig — Lmin/S —0.33 fm (0.5 fm) in the ground 
state of -B. As noted already, the energy increase due to the string stretching 
is of the order of 1 GeV/fm. The string energy in B is Sruq ~ crLrain — 1 GeV 
(1.5 GeV), which leads to mj ~ (—0.5 GeV). The color Coulomb force gives 
the moderate attractive energy of about -0.2 GeV (-0.1 GeV).^) The effect 
is much weaker than those (~ —0.5 GcV) obtained in calculations wiht use of 
the quark cluster model and the bag model, because the probability of close 
approach of the quarks is much smaller for the Y-shaped configuration. Energy 
contribution to mj, being ~ 0.2 GeV (—0.4 GcV), is to be provided as energy 
concentration at the junction. Using the relation 2mq ^ mij, we can estimate 
mij ~ 0.7 GeV (1 GeV), the energy needed for formation of an interjuction. 
Thus there remains a problem how to understand the value of mj dependent on 
the choice of Lmin- As mentioned in §3, following the lattice QCD calculations, 
the case of Lmin — 1-5 fm seems favorable, since for Lmin — 1-0 fm the string length 
£s is comparable with the width of the flux tube and the string-junction picture 
looks not so apparent. For the case of Lmin — 1-5 fm, the negative energy of the 
order of mj ~ —0.4 GeV is inevitable, since the string energy is considerably larger 
than niB, but this feature is not unreasonable in the following sense. The negative 
mj means the concentration of negative energy at the junction in the ground state 
of B without interj unction. In the exotic baryon states, however, the formation of 
new junctions inevitably accompanies the creation of interj unction (s) which needs 
the large positive energy in unit of m/j ~ 1 GeV. 

Although the value of mj looks directly linked to the parameter Csq in Eq. (2-2), 
estimate of mj from C^q is not possible because C-^q contains an integration constant. 
Various effects to the energy of the ground state of B not specified in the simple mass 
formula in the SJM would contribute to mj, and this problem is beyond the scope 
of the present study. As noted already this problematic point is not obstacle in the 
present work, because we have used the quantity S preferentially. 

6.2. On the BB repulsive core 

The NN repulsive core is one of the main factors leading to the property that 
the nucleus with the low density is the lowest state in hadronic matter. On the BB 
repulsive core, there have been rarely studied by taking into account the confinement 
mechanism, which is the nonperturbative effect originating from the color degrees 
of freedom in QCD. In the studies in the quark cluster model, although the quark 

In the 1980's there was pointed out the possibiUty that the strange quark matter with almost 
equal numbers of u, d and s is the absolutely lowest state of hadronic matter. In a previous study 
with use of the quark cluster model and the bag model,""^' we have shown that this possibility is 
not realized. 
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confinement is incorporated in a potentilal form between quarks, the interaction 
between clusters (baryons) is described in such a way that the confinement plays no 
substantial role as far as the color singlet is imposed to the total system, in spite of 
its significant contribution to the internal energy of the cluster of quarks (baryon). 
Since these studies have been developed to a stage that most of the existing data of 
the BB system are well explained,^^)'^^^ the important findings are surely contained 
in the results of these works. Neverthless it seems strange that the effects of the 
confinement docs not appear in any essential manner in the BB interaction. We feel 
that something essential may be lacked in the recognition hihterto obtained. In the 
present paper we have shown a straightfoward way to understand the BB repulsive 
core as originationg from the confinement. 

Recently the first realistic approach on nuclear force in lattice QCD calculations 
has been made in the quenched approximation.'^^-' Although the outside attraction 
and the inside repulsion of the NN S-wave interaction have been reported, it is 
premature to discuss the results in quantativc level, and wc expect future progress. 
From the viewpoint of the present paper, it is desirable to see at what extent the 
process like JJ-pair creation takes part on the short-range interaction, although it 
will be a task beyond the quenched approximation. This problem is connected with 
the study to clarify whether exotic dibaryons exist or not. It seems quite natural 
that the lattice QCD simulation has not been performed yet concerning the exotic 
dibaryons, since they appear through the virtual creation of the BB pair and that 
of the qq pairs, following the context of the present paper. To get observational 
evidence on these exotic dibaryon states, we expect future experiments. Concerning 
the dibaryon problems, it seems promising to obtain exclusive data with use of the 
anti-proton beam such as p -|- ^He — > + 7 and p + ^He 0^ + N. 

6.3. On three-baryon interaction 

Hitherto in most cases, the three-body interaction in the BBB system has been 
introduced as a supplementary effect to reproduce the data in a quantitative level. 

In recent years, its indispensable role in the nuclear system has been recognized 
progressively, concerning the saturation problem, the maximum mass of neutron 
stars, and the requirement to reproduce various kinds of precise measurements on 
few-miclcon systems. Nevertheless, it seems fair to ask why thrcc-baryon interaction 
plays a decisive role in the fundamental aspects of the nucleus. We can accept that 
its importance is natural, if its short-range part originates from the confinement and 
its intermediate part comes from meson exchange, especailly pion exchange. 

6.4. On exotic multi-baryon system 

In this paper we have treated the exotic baryon systems of the baryon num- 
bar N = 2, 3 in the SJM. What aspects of such exotic systems come about when 
N becomes much larger? We can apply the same procedure as shown in Fig. 6 to 



*' Recently Imachi, Otsuki and Toyoda have studied unconventional hadrons including exotic 

ones, extending the string-junction model they developed in the 1970's.'"^-' Further study by assign- 
ing the recently observed peak at the mass 4.43 GeV {Z^)^^> as the exotic meson M| (e.g., udcc for 
Z+), has suggested the appearance of the dibaryon Dg here discussed at the mass about 4.3 GeV.*^^ 
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the N = A system, and find that its energy is higher than the ordinary value 4ms 
by ~ 3f . Furthermore we can extend the treatment to the general case. Such an 
exotic multi-baryon system (abbreviated to EMBS), whose possibility was pointed 
out in Rcf. 7c), have the energy higher than NniB of the ordinary baryon system 
by (N — 1)£. Namely, EMBS has the saturation property of energy; £^/baryon~ 2 
GeV/baryon. However to estimate the saturation density of EMBS is problematic, 
contrary to the unifom matter. The reason is that EMBS connected by strings, junc- 
tions and interj unction strings is flexible in spatial form like the "chain". Therefore 
it seems difficult to consider EMBS beyond few-baryon systems, unless its spatial 
configuration is well defined. 

§7. Concluding remarks 

In this study we have shown a way to unifiedly understand the universal repulsive 
core of the BB interaction and the universal repulsion of the three-baryon interac- 
tion, standing on the the string-junction model where the confinement mechanism 
in QCD is built in. 

The key points of recognition are summarized as follows. In order for the con- 
finement mechanism to persisit at closest approach of baryons, the BB-paiv (more 
widely the pair of the junction and anti-junction) excitation takes place and leads 
to the formation of a new closed string-junction net with exotic nature. For the BB 
system the energy needed for the BB-paii creation appears as the rcpTilsivc core 
with ~ 2 GcV height. For the BBB system a half of the energy for the 2BB-pa,iT, 
being ~ 2 GeV, is regarded as the strength of the three-body potential. Partici- 
pation of such exotic baryonic states brings about the latent effect on the ordinary 
baryon states, which appears as the short-range repulsion. Thus we can understand 
the universal (flavor- and spin- independent) short-range repulsion as originating 
from the confinement mechanism. Results of numerical calculations show that the 
universal repulsion of three-baryon interaction derived in the string-junction model 
provides us with the sufficiently stiff equation of state, which is consistent with the 
observations of neutron star masses, even when the hyperon-mixing occurs in the 
neutron star core. 

We have developed our studies on the assumption of the existence of exotic 
multi-baryon states. Hitherto experimental evidence about exotic hadrons has not 
been obtained yet. Since there is no reason in QCD to deny the existence of hadrons 
other than the ordinary ones, this is a basic problem to be solved in future in hadron 
physics. Recently exotic hadrons have attracted much interest, exprimetally and 
theoretically,^^) and we hope that active studies now being in progress will elucidate 
such basic problems in future. 

Concerning the problem under investigation, the reasoning taken here will not 
lose validity even though exotic hadron states with narrow width are not observed. 
The reason is that, e.g. for the BB case, if there exist a group of states with a 
junction number differnt from that of the ordinary BB case almost restricted to the 
small spatial space corresponding the core region of BB and their energies are as 
high as 2 GeV, the ordinary BB system can scarcely enter into such region. This 
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implies the existence of the repulsive core. 
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